Phonon-assisted tunneling in asymmetric resonant tunneling structures 
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Based on the dielectric continuum model, we calculated the phonon assisted tunneling (PAT) cur- 
rent of general double barrier resoant tunneling structures (DBRTSs) including both symmetric and 
asymmetric ones. The results indicate that the four higher frequency interface phonon modes (es- 
pecially the one which peaks at either interface of the emitter barrier) dominate the PAT processes, 
which increase the valley current and decrease the PVR of DBRTSs. We show that an asymmetric 
structure can lead to improved performance. 
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I. INTRODUCTION 

The double barrier resonant tunneling structure 
(DBRTS) H continues to attract attention both for its 
potential applications to electronic devices and for its 
value in exploring fundamental phenomena p|, including 
tests of electron-phonon coupling theories [p|-]l3[ . Despite 
extensive study of the DBRTS, challenging problems con- 
tinue to exist with respect to practical applications and 
also for modeling the characteristics of the DBRTS . 
Our objective here is to fully analyze the phonon-assisted 
tunneling (PAT) current, especially in the valley region, 
and to demonstrate that there can be important advan- 
tages in designing a DBRTS with an asymmetric struc- 
ture. We include all of the phonon modes and provide 
a complete analysis of the relative importance of these 
modes, particularly for the asymmetric DBRTSs where a 
limited investigation has been performed to date |l2|,|l3| ■ 
It is important to consider all of the phonon modes when 
the asymmetric DBRTS is considered 14 16| . 

The motivations for considering the asymmetric 
DBRTS are the possibility for improving performance, 
tunability of the negative differential resistance (NDR) 
p7|, and tuning of charge accumulation in the quantum 
well (QW) [|l8|,^9). This asymmetry can be introduced 
by varying the material compositions in the two barriers 
or by creating two barriers of unequal width, and we will 
consider both possibilities here. Our treatment includes 
a detailed study of all of the phonon modes, assuming the 
dielectric continuum theory and incorporating the effects 
of subband nonparabolicity. 

Our main concern is with DBRTS performance. A 
large NDR is particularly important for high-frequency 
resonant tunneling and rapid switching devices and is 
usually quantified by the peak-to-valley ratio (PVR) of 
the current-to- voltage characteristic curve. A large PVR 
and low valley current are desirable for most resonant 
tunneling diodes (RTD) device applications. Generally, 
PAT, r — X intervalley tunneling, impurity scattering, 
the interface roughness scattering and the tunneling of 



quasi-two-dimensional subband electron in the pseudo- 
triangular well in the emitter can cause the valley cur- 
rent [[10|. For a polar semiconductor DBRTS, the effects 
of phonon scattering on resonant tunneling are very im- 
portant and inevitable especially at room temperature. 
The electrons in the DBRTS may emit phonons during 
the resonant tunneling process. In general, the contribu- 
tion by the PAT current is small compared with the co- 
herent resonant tunneling process but can be large com- 
pared to the other contributions to current. Hence the 
total tunneling current density can be approximated by 
J = J c + J p , for J c the coherent tunneling current den- 
sity and J p the PAT current density. Although the PAT 
current has been investigated theoretically p- Q , [l3| , we 
provide a more rigorous investigation of electron-phonon 
scattering and the PAT current in general DBRTS in- 
cluding both symmetric and asymmetric ones, with the 
contribution of all of the phonon modes accounted for 
quantitatively. Hence the PAT physical picture is fur- 
ther clarified. 



II. THEORY 

Working within the framework of the dielectric con- 
tinuum model and treating the electron-phonon inter- 
action via the Frdhlich-like Hamiltonian, we can calcu- 
late the electron-phonon scattering rate W according to 
the Fermi golden rule for two cases jL5|: scattering by 
interface phonons W-^*^' (fcj, E z ) and by confined LO 

phonons Wj^*^ (fcj, E z ). The scattering rate depends 
explicitly on the phonon occupation number iV p h which 
is temperature-dependent according to the Planck dis- 
tribution p2| . As the width of the final resonant state 
is very narrow, the final state is treated as a com- 
pletely localized state in the well ||-|,|Jl|Q. We have 
recently presented the expressions of these scattering 
rates ^ntLo(ki) E z ) for both symmetric and asymmetric 
DBRTSs 01. 
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Electron tunneling in a DBRTS depends sensitively on 
the bias voltage V. The coherent tunneling current den- 
sity J c and PAT current density J p are given by 



and 



J c = en(T(E z )v z (E z )) 



J p = eN(W)/A. 



(2.1) 



(2.2) 



Here e is the absolute value of the electron charge, N 
the total electron number in the electron reservoir with 
volume f2, n = N/Cl the electron density (assumed to 
be constant), T{E Z ) the electron transmission coefficient, 
and v z (E z ) the longitudinal electron velocity, A the cross- 
sectional area of the structure, and W the electron- 
phonon scattering rate fl2]| . () represents averaging on 
quantum stat es. 

From Eq. (2.1) we obtain M 



J c Jr — > Jr. 



(2.3) 



with 



Jr.^ = 



T(E Z ) 



27r 2 r jo 

ln{l + exp[(E F - E z )/k B T]}dE z 



(2.4) 



and 

Jc- 



/ T(E Z ) 
Jo 



em\\kBT 

ln{l + exp[(£ F - eV - E z )/k B T]}dE z 



(2.5) 



Here J c ^ {J c ^) denotes the tunneling current density 
from the emitter (collector) to the collector (emitter), 
and J c ^ is critical for the NDR of the DBRTS, which 
controls the coherent tunneling current density. Ep and 
Ep — eV are, respectively, the local Fermi energy levels in 
the emitter and collector, and my is the electron effective 
mass in the x — y plane parallel to the interfaces of the 
DBRTS. 



From Eq. (2.2) we derive 
emukpiTL 



Jp — 



1/2 



E 
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x In [l + e (E F -E z )/k B T \ Wp ( Ez } dE: 



e: 



-1/2 



(2.6) 



with L e the emitter length, m z the electron effective mass 
for longitudinal motion in the emitter, v — 2, 3, 4 the in- 
dex for confined bulk-like LO phonon modes in the left- 
barrier, the well and the right-barrier layers, respectively, 
and m (1 to 8) the index for the eight interface phonon 
modes. The electron-phonon scattering rate in the v th 
layer is Wv(E z ), and W m (E z ) is the scattering rate be- 
tween electron and the m th interface phonon mode. Ex- 
plicit expresssions for W„(E Z ) and W m (E z ) have been 
presented in our recent paper [12]. 



In obtaining Eq. (2.6) we have assumed W(ki,E z ) = 
W(0, E z ) = W(E Z ) and have ignored the contribution 
from collector-to-emitter tunneling via a phonon absorp- 
tion||f]. 

De spit e the apparent minor difference between 



Eqs. (|2.4| ) and ([2J]), simulating J p is more costly in com- 
puter time compared to simulating J c because W P (E Z ) 
is more complicated than TiE z ) [[t2|. Our calculations 
show that J p mainly comes from the electron scattering 
by the higher frequency interface phonons (especially the 
interface phonons localized at either interface of the left 
barrier). 



III. NUMERICAL RESULTS AND DISCUSSION 

Numerical calculations have been performed for asym- 
metric and symmetric DBRTSs A{x, d) defined as 

A(x,d) = u + GaAs(1000l)/Al x Ga 1 _ a As(30l)/GaAs(60A) 



/Alo. 3 Gao. 7 As(d)/n + GaAs(1000A), 



(3.1) 



with my = m z and equal doping concentration 10 cm 
in emitter and collector. The physical parameters used 
are the same as in |l(| . 

Figure |l| p resents the dispersion calculated from 
Eq. (4) of ]15| ] for the eight interface modes of structure 
,4(0.25, 204.), which reveals that the four lower-frequency 
modes occupy a much narrower frequency band than the 
four higher-frequency modes. Moreover, the dispersion of 
the interface modes is significant for the case k < O.lA -1 
and negligible for k > O.lA^ 1 . The electron-interface- 
phonon coupling function T(k, z) [ [l5[ is a complicated 
function of both z and k in a DBRTS. The T(k, z) - z 
relation in Fig. |2| reveals the strength of the electron in- 
teraction with different interface modes peaks at differ- 
ent interfaces. For example, the electron interaction with 
mode 7 (denoted e-p(7)) peaks at the z = and 30 A in- 
terfaces (i.e., either interface of the left barrier). We also 
find that |r(fc, z)\ decreases rapidly as a function of k for 
< k < O.OlA" 1 , slowly for O.Oll^ 1 < k < 0.081" 1 , 
and then rapidly for k > 0.08A -1 for the e-p(7) inter- 
action. The 7 th interface mode is much more important 
than the other modes in the asymmetric DBRTS, and 
the four higher-frequency modes produce intensive po- 
larisation in the DBRTS, resulting in a significant inter- 
action with electrons. On the contrary, the four lower- 
frequency modes produce a weak interaction with elec- 
trons compared with the higher-frequency modes, which 
can be ignored [^2| , p4| -[l6| . In the following calculations, 
we will thus only consider the contribution of the four 
higher-frequency modes to the scattering rate and the 
PAT current for simplicity. 

Figure || shows the scattering rate W vs incident elec- 
tron energy E z for ,4(0.25, 20l). For n = 10 18 cm" 3 , the 
Fermi energy level E F = 42.5 meV at T = 300 K. We 
observe that the contribution due to interface phonons 
is larger than that due to LO bulk-like phonons for the 
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electron with lower incident energies. We know from the 
Fermi distribution function that the emitter states are 
appreciably populated only for E z < Ep + ksT = 68.3 
meV (as T = 300 K). Hence, we can infer that the in- 
terface phonons contribute much more than the confined 
bulk- like LO phonons to the PAT current. 

We present the electron-phonon scattering rate, in- 
clouding the subband nonparabolicity, as the dash-dot- 
dot line in Fig. [| which shows that the subband non- 
parabolicity has a large influence on the electron-phonon 
scattering. The peak position shifts to a lower energy, 
and its value decreases under the influence of the sub- 
band nonparabolicity. 

In the case of including subband nonparabolicity, PAT 
current-to- voltage curves are shown in Fig. |i| at T = 300 
K for the structure ,4(0.25, 20A). Figure 0(a) shows the 
tunneling current assisted by the four higher-frequency 
interface phonon modes and their sum. We can see from 
Fig. |](a) that the 7 th interface mode, which is localized 
at either interface of the left barrier, is the most im- 
portant of all of the interface modes, and this result is 
consistent with the results shown in Fig. || The total 
interface PAT current is a complicated function of the 
applied voltage and has two peaks for increasing voltage. 
Figure ||(b) gives the confined bulk-like LO PAT current 
density in structure .4.(0.25, 20A). This figure indicates 
that the PAT current from the LO phonons in the well is 
much larger than those from the LO phonons in the two 
barrier layers with a complicated behavior for increasing 
bias voltage. Figure |J(c) presents the total PAT current 
density including the interface and the confined bulk-like 
LO phonons. We can see from Fig. [|(c) that the inter- 
face PAT current is larger by one order of magnitude than 
the confined LO PAT current, confirming that interface- 
phonon scattering dominates over confined LO phonon 
scattering (c.f. Fig. 3). Moreover, Fig. 4(c) also shows 
that the total PAT current is a very complicated function 
of the applied voltage and has two peaks, similar to the 
results based on the Green function method This is 
completely due to the complexity of the contribution of 
the phonon modes to the PAT current in our DBRTS, as 
shown in Figure 4(a). Figure || also indicates that the 
PAT current is mainly determined by scattering between 
electrons and higher frequency interface phonons (espe- 
cially the interface phonons localized at either interface 
of the left barrier), showing a clear physical picture for 
the PAT process in general DBRTS. 

Figure g shows the total current-to-voltage curve at 
room temperature for .4(0.25,204), including coherent 
and PAT currents. This figure shows that PAT increases 
the valley current and decreases the PVR. The result 
shown in Fig. 5 is similar to those obtained by the Wigner 
function method Hi and the Wannier function envelope 
equation method fllpf . 

We have also studied the current-voltage characteris- 
tics for DBRTSs .4(0.3, d) with d = 20, 30 and 40A. The 
calculated results show that when the right-barrier thick- 
ness d increases, the peak current decreases, the PVR 



increases, and the peak position shifts towards higher 
bias voltage. These three characteristics are in agree- 
ment with recent experimental results Jl8| |. Moreover, 
we have also studied the current-voltage characteristics 
for A(x, 30A) with x = 0.2, 0.3 and 0.4. The calculated 
results show that when the Al composition x, and thus 
the height of the left barrier, increases, the peak current 
decreases, the PVR increases, and the peak position re- 
mains at the same bias voltage. These theoretical results, 
which show that an asymmetric DBRTS can lead to im- 
proved performance, await experimental confirmation. 

IV. SUMMARY 

Employing the dielectric continuum model, including 
all the phonon modes and treating conduction band non- 
parabolicity, the PAT process is investigated in detail. 
We show that for a DBRTS, no matter it is symmet- 
ric or asymmetric one, the four higher frequency modes, 
which reduce to two symmetric and two antisymmetric 
higher frequency modes for a symmetric structure, dom- 
inate the interface PAT process. In particular, for the 
interface PAT process in a symmetric DBRTS, the two 
symmetric higher frequency modes are most important, 
the two antisymmetric higher frequency modes are less 
important, and the two symmetric and antisymmetric 
lower frequency modes are negligible. The above opinion 
is different from that of Refs. 0-0]. 

In general, the confined LO phonons in the well layer 
are more important than those in the two barrier layers, 
and the four higher frequency interface modes (i.e., mode 
5 to 8, especially the 7 th one which peaks at either inter- 
face of the left barrier) dominates over the four lower 
frequency interface modes and all of the confined LO 
phonon modes to electron-phonon scattering and PAT 
current. The PAT increases the valley current and de- 
creases the PVR of DBRTS. It is worth mentioning that 
the mode 7 becomes the highest frequency symmetric in- 
terface mode for a symmetric DBRTS. The PAT physical 
picture stated in the above is useful and important for 
further understanding PAT process in DBRTSs and for 
designing better RTD devices. 

Subband nonparabolicity has a significant influence 
on electron-phonon scattering, PAT, and the current-to- 
voltage characteristic of a DBRTS. 

We find that the peak current is reduced, the position 
of peak current is shifted to a higher voltage, and the 
PVR is enlarged if the right-barrier width is increased 
when the two barriers have the same height. The peak 
current is reduced and the PVR is increased by suitably 
increasing the left-barrier height when the two barriers 
have the same width. An asymmetric DBRTS with a 
suitably designed structure has a larger PVR than its 
commonly- used symmetric counterpart. The results ob- 
tained in this paper are useful for analysing coherence- 
breaking phonon scattering and for potentially important 
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resonant tunneling device applications. 
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FIG. 2. Spatial dependence of the normalized coupling 
functions (he 2 / Ae y 1/2 V(k, z) for the structure .4(0.25, 2uA) 
(k — 0.014 -1 ). Here the numbers by the curves represent the 
interface-phonon frequency in order of increasing magnitude: 
(a) for the four lower-frequency modes and (b) for the four 
higher- frequency modes with the interfaces localized at z =0, 
30, 90 and 110 A, respectively. 

FIG. 3. The normalized electron-phonon scattering 
rate W/(N p h + 1) vs incident electron energy E z for 
.4(0.25,201) at the bias 150 mV. The dashed line and 
dash-dot line represent, respectively, the contribution of the 
interface phonons and confined bulk-like LO phonons, and the 
solid line is their sum in the absence of subband nonparabolic- 
ity. The dash-dot-dot line is the total scattering rate including 
subband nonparabolicity. 
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FIG. 4. PAT current-to-voltage curves at room tempera- 
ture for .4(0.25, 204), including the subband nonparabolicity: 
(a) Interface PAT: the dashed line represents the 5 th interface 
mode contribution for the tunneling current, the dash-dot-dot 
line for the 6 th mode, the thin solid line for the 7 th mode, 
the dash-dot line for the 8 th mode and the heavy solid line 
is their sum. (b) LO PAT: the dashed line represents the 
confined LO phonon in the left barrier (Alo.25Gao.75 As) con- 
tribution for the tunneling current, the dash-dot line for the 
LO phonon in the right barrier (Alo.3Gao.7As), the thin solid 
line for the confined bulk-like LO phonon in the well (GaAs) 
and the heavy solid line is their sum. (c) Interface PAT cur- 
rent (thin solid line), confined LO PAT current (dashed line) 
and their sum (heavy solid line). 

FIG. 5. The total current-to-voltage characteristic curve 
at room temperature considering the subband nonparabolic- 
ity for .4(0.25, 204). The dashed line stands for the coherent 
tunneling current density. The solid line is the total tunnel- 
ing current density combining coherent tunneling and PAT 
currents. 



FIG. 1. The dispersion curves of the interface modes for 
structure .4(0.25, 20A). 
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